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Abstract 

An appropriately oriented D3-D7-brane system is the holographic dual of relativistic 
Fermions occupying a 2+1-dimensional defect embedded in 3+1-dimensional spacetime. 
The Fermions interact via fields of Af = 4 Yang-Mills theory in the 3+1-dimensional 
bulk. Recently, using internal flux to stabilize the system in the probe A7 << A3 limit, 
a number of solutions which are dual to conformal field theories with Fermion content 
have been found. We use holographic techniques to study perturbations of a particular 
one of the conformal field theories by relevant operators. Generally, the response of a 
conformal field theory to such a perturbation grows and becomes nonperturbative at 
low energy scales. We shall find that a perturbation which switches on a background 
magnetic field B and Fermion mass m induces a renormalization group flow that can be 
studied perturbatively in the limit of small m 2 / B. We solve the leading order explicitly. 
JL . We find that, for one particular value of internal flux, the system exhibits magnetic 

qj catalysis, the spontaneous breaking of chiral symmetry enhanced by the presence of the 

magnetic field. In the process, we derive formulae predicting the Debye screening length 
of the Fermion-antiFermion plasma at finite density and the diamagnetic moment of the 
ground state of the Fermion system in the presence of a magnetic field. 



1 Introduction 

The AdS/CFT correspondence pQ offers the hope of direct, mathematically precise 
and systematically correctable study of the strong coupling limit of some quantum 
systems [2]. Condensed matter physics in particular encounters a number of systems 
which exhibit quantum critical behavior and where the coupling can be argued to be 
strong. In this Paper, we shall study the holographic dual of 2+1 dimensional quantum 
field theories with relativistic fermions. Potential applications could be to condensed 
matter systems which have emergent relativistic 2+1-dimensional Fermions, examples 
of which are graphene [3] [3], topological insulators [3], the D-wave state of high T c 
superconductors [3] and simulation of such systems on optical lattices [7] . 

The Coulomb force in graphene in particular is strong. However, it also violates 
the relativistic Lorentz symmetry of the free low energy electrons. Our analysis in the 
following, being relativistic, only applies if graphene finds a way to be relativistic even 
in the presence of strong nonrelativistic interactions. This could happen, for example, 
if the relativistic theory is a conformal field theory occurring at the infrared fixed point 
of a renormalization group flow. There are some experimental indications that this 
could be the case. However, there is little theoretical support for this idea, part of the 
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difficulty being the absence of reliable techniques for the strong coupling regime. We 
will not address this problem directly in this paper. What has been done in previous 
work [5] is to demonstrate that, at strong coupling, conformal field theories that are 
viable candidates do indeed exist. Here we shall examine some of the properties of the 
strongly coupled conformal field theory. In particular, we shall be interested in the fate 
of the field theory when it is perturbed by relevant operators. Such perturbations, such 
as turning on finite charge density, external magnetic fields or a Fermion mass operator 
corresponding to sublattice asymmetric charge density are very relevant to the physical 
properties of such systems. 

In weakly coupled field theory, the propensity of gauge field mediated interactions to 
form a chiral condensate and break chiral symmetry is greatly enhanced by the presence 
of an external magnetic field. This phenomenon is called magnetic catalysis [5]-|17|. 
The interesting question as to whether it persists at strong coupling has been addressed 
in some holographic models where it has indeed been found to occur, particularly in 
holographic D3-D5 systems [H]-|123- However, but it has proven to be more elusive in 
the D3-D7 system [5] and there are even cases of "anti-catalysis", suppression of a chiral 
condensate by a magnetic field [2HIEI]- Here, with our explicit perturbative solution 
of the holographic system, we shall find that magnetic catalysis can indeed occur, but 
only for one special value of a particular tuneable parameter. 

An example of a top-down holographic construction of strongly interacting 2+1- 
dimensional relativistic Fermions uses appropriately oriented probe D7-branes in the 
AdS 5 x S 5 geometry that is sourced by N coincident D3-branes [32] QU [33] . The 
holographic construction begins with D7 and D3 branes oriented as in Table 1 (jl.ll) . 
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Table 1 : D3-D7 orientation 



The N Z)3-branes and the -/V7 D7 branes are extended in 2+1-spacetime dimensions 

2; 1 , x 2 ) = (t,x,y) where they have 50(2, 1) Lorentz symmetry. The lowest energy 

states of the 3-7 open strings are N species of 2+1-dimensional 2-component Fermions. 

The x 9 direction is orthogonal to both the D3 and D7. The D3 and D7 can be separated 

in that direction, introducing a bare mass for 3-7 strings. For 2-component Fermions, a 

bare mass must violate parity. We will discuss how parity is formulated in the D-brane 

construction shortly and we will see that parity must be formulated to change the sign 

of the separation of the D3 and D7 branes. 

To apply holography, the limit where the number of jD3-branes N is large is taken 

while holding the product of N and the closed string coupling constant g s fixed. In 

the holographic duality, the closed string coupling constant is related to the Yang-Mills 

coupling of the bulk AT = 4 supersymmetric Yang-Mills theory by 4:irg s = g YM - The 

quantity which is held fixed in the large N limit is A-Kg s N = gy M N = A, the 't Hooft 

coupling of the gauge theory. Then, the D3-branes are replaced by the AdS& x S 5 

1 

geometry. The radii of curvature of the AdS§ and S 5 are L = A 4 Vff . The Z?7-branes 
are treated as probes and the dynamical problem is to find their embedding in AdS§ x S 5 . 

This D3-D7 configuration has a unique feature that it is non-supersymmetric, but is 
free of tachyons and the only low energy modes of the D3-D7 open strings are Fermions. 
As a consequence, the decoupling limit produces a system which at weak coupling 
contains only chiral Fermions. Being a non-supersymmetric configuration, the D3 and 
D7-branes repel each other. This shows up as an instability that appears when one 
attempts to embed the D7-brane in AdSr> x S 5 . Fluctuations of the geometry violate 
the Breitenholder-Freedman bound in the large AdS radius regime. This instability can 
be fixed by introducing flux of the world-volume gauge fields of the D7-brane, either 
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Figure 1: Defect conformal field theory: Fermions are constrained to occupy a plane, denoted by the 
vertical line through the center of the diagram. This plane divides the three dimensional space into two 
regions which are occupied by four dimensional conformal M = 4 supersymmetric Yang-Mills theories with 
different gauge groups. 

an instanton bundle [TH] or U(l) magnetic magnetic monopole fluxes [33]. In the latter 
case, which is the one we will focus on in this paper, the four dimensions of the D7 
world- volume which are embedded in S 5 are taken as two 2-spheres, S 2 and S 2 , and 
each 2-sphere has a number njj and no units of Dirac magnetic monopole flux. It 
was shown in reference [33] that the latter configuration is stable, at least to small 
fluctuations, if either njj or fiu is large enough. 

The quantum field theory which is dual to the D3-D7 system is a defect field theory 
consisting of Fermions confined a 2-dimensional plane which separates 3-dimensional 
space into two regions as depicted in figure [T] There are Nj species of 2-component 
spinors of the 2+1-dimensional Lorentz group 50(2,1). In the conformal invariant 
solution this symmetry is extended to the 2+1-dimensional conformal group 50(3,2). 
The Fermions transform in the fundamental representation of a global U (N7) symmetry. 
Since we consider no processes which use the non-abelian nature of U(N^), for simplicity 
we will take N7 = 1 (and remember that to apply to graphene, we need Nj = 4 to 
produce the correct flavor symmetry). The Fermions also transform in the fundamental 
representation of the gauge group of Af = 4 supersymmetric Yang-Mills theory which 
inhabits the 3+1-dimensional bulk of the space-time. The bulk Yang-Mills theory has 
different gauge groups on each side of the defect, as shown in figure [1] This is a result of 
the fact that, the D7-branes with internal fluxes which we shall use can be described as 
D7-branes with noun D3-branes dissolved into their worldvolumes. The D7-brane then 
forms a boundary between regions with different numbers of D3-branes and therefore 
different amounts of Ramond-Ramond 4-form flux, thus different ranks of the gauge 
group in the field theory dual. 

We shall be interested in field theories which become parity and charge conjugation 
invariant in their high energy limit. This is what is expected in a class of condensed 
matter theories where the underlying dynamics is parity and particle-hole symmetric. 
These symmetries can be broken by deformations like external magnetic field, chemical 
potential or parity violating mass terms which are irrelevant in the ultraviolet limit 
but can be important to the infrared properties of the theory. On the D-brane side, 
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imposing parity and charge conjugation symmetry will involve taking the appropriate 
boundary condition for the embedding of the D7-brane in AdSs x S 5 as well as setting 
the fluxes equal, rip = fin- 

The paper is organized as follows: In Section 2 we review the embedding of the 
probe D7-brane in AdSs x S 5 . In Section 3 we discuss the properties of the conformally 
invariant solution of the embedding problem. In Section 4 we discuss the solution with a 
chemical potential and charge density. In particular, we derive expressions for the Debye 
screening length at strong coupling, as functions of chemical potential and of density. In 
Section 5 we examine the same system with an external magnetic field in the special case 
that the charge density is tuned to zero. We find a simple expression for the diamagnetic 
moment of the system. We solve the embedding equation perturbatively in the ratio of 
condensate to magnetic field. We find a relationship between the mass m and the chiral 
condensate c to linear order in m 2 /B in equation (|5.96p . In Section 7 we show that 
turning on an infinitesimal charge density can also be taken into account perturbatively 
and we write the embedding equation to leading order in the filling fraction p/B. 



2 D7-brane 

In the limit where the string theory is classical, the problem of embedding a D7-brane in 
the AdS$ x S 5 geometry reduces to that of finding an extremum of the Dirac-Born-Infeld 
and Wess-Zumino actions, 
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- y/- dct(g + 2na'F) + C< 4 > A F A F 



(2.2) 



where g s is the closed string coupling constant, which is related to the Af = 4 Yang-Mills 
coupling by 4:irg s = g\ u -, {(Jo, <xi, . . . , 177} are the coordinates of the D7-brane world- 
volume, g a b{<?) is the induced metric, is the 4- form of the AdS$ x S 5 background, 
F is the world- volume gauge field and the D7 brane tension is 



T 7 
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(2.3) 



(27r) 7 a' 4 

We shall work with coordinates where the metric of the background space AdS§ x S 5 is 



ds 2 



L 



\-dt 2 +dx 2 + dy 2 
+dip 2 + cos 2 ip(d8 2 + sin 2 



, 2n dr 2 
dx 2 ) + -^r- 



sin 2 iP{d9 2 



sin 



(2.4) 



Here, (t, x, y, z, r) are coordinates of the Poincare patch of AdS$. In our notation, and 
in natural units h = 1 and c — 1, r has the dimension of inverse length and (x, y, z,t) 
have dimensions of length. L is the radius of curvature. The boundary of AdS§ is 
located at r — > 00 and the Poincare horizon at r — > 0. The 5-sphere is represented by 
two unit 2-spheres, S 2 with polar coordinates (8, 4>) and S 2 with coordinates (8, <fi). The 
2-spheres are fibered over the interval ip G [0, V\. The Ramond-Ramond 4-form is 



C (i) = £ Vdi A dx AdyAdz + L 



d cos 8 A dd> A d cos 8 A > 



with c(i/)) obeying the equation 



8 sin 2 -0 cos 2 4> 



(2.5) 



(2.6) 



The dynamical variables are the ten functions of eight world- volume coordinates which 
embed the D7-brane in AdS 5 x S 5 , 



{x(a), y(a),z(a),t(a), r(a), j>(a), 8{a), 0(a) , 8 \a) , 4>(a)} 
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as well as the eight worldvolume gauge fields 

{A (a),A 1 (a),...,A 7 (a)} 

Parity and charge conjugation are important symmetries for the class of quantum 
field theory systems that we are interested in. For example, graphene is certainly parity 
invariant and, to a good approximation, it has particle-hole symmetry. We would expect 
to model it using field theories that have parity and charge conjugation symmetry. The 
string theory dual of such a field theory should also have these symmetries. We should 
therefore make sure that the problem of finding a minimum of the action (|2.2p itself is 
symmetric. Parity in two dimensions is a reflection of one of the spatial coordinates, 
x — > — x. This symmetry is usually broken by Wess-Zumino terms, which ar(£] 

j d & ae^-^d^t(a)d^x(a)d^y(a)d^z(a)r A (a) d^A^(a)d^ 7 (a) (2.7) 

J (fae^-^O^ cos 0(0")^ 0(a) cos9(a)d^(a)d (^^A^d^A^a) 

(2.8) 

Indeed, (|2.7[) changes sign when x — > —x. We must therefore compensate the sign 
change by another change of the world-sheet variables, say o\ — > —<T\. However, now 
(|2.7p is invariant but (|2.8p changes sign and is not invariant. We can also make it 
invariant by another change of variables, (ip, 0, 0, 0, 0) — > (S- — ip, 0, 0, 6, 0). (Note that 
c'(ip) is invariant and dip changes sign under this transformation.) Then, both of the 
Wess-Zumino terms are invariant, as is the Dirac-Born-Infeld action. In summary, our 
parity transformation is the replacement 

= {-x(a), y(a), z(a), t(a), r(a), f - ^(a),§(a), 0», 0(a), 0(a)} (2.9) 
4>') = -AiW, A 2 (a), A 7 (a)} (2.10) 

^' = {a , -cri,^,..., <7 7 } (2.11) 

Charge conjugation symmetry (C) is the replacement A^ — > — A^ and the Wess-Zumino 
terms are invariant. However, we shall introduce a background field 

F = (dcos0 A d(f> + d cos A d(j>) (2.12) 

27ra' 2 V / 

Here, / is the strength of the monopole bundl^l which is needed to stabilize the system. 
Note that, to be invariant under parity, (|2.9p - (|2.1ip . the fluxes on the 2-spheres have 
to be equal, nu — no- This is seen most clearly by noting that, in (|2.9[) - (|2.1ip , the 
transformation interchanges the 2-spheres. The background field breaks a symmetry 
if the two components in the total field F = Fq + 5F transform differently. Under 
C, 5F —> —8F. We shall need a definition of C so that, that F — s- — F . Such a 
transformation is 

C : {x'(a),y'(a),z'(a),l/(a)y(a),ij'(a) > e'(a) > cl>'(a),e'(a),^(a)} = 

= {x(a), y(a),z(a),t(a), r(a), i&(tr), 9(v),2* - 0(a), 0(a), 2tt - 0(a)} (2.13) 

4» = {-M<r),-M<r),-M<r), • • • . ~M<r)} ( 2 - U ) 

x Note that we have presented the second Wess-Zumino term (|2.8p in a form that is integrated by parts. 
This was to avoid specifying the integration constant in c(ip) which would be obtained by integrating the 
expression in (|2.6[) . Specifying that integration constant (which we shall do in the following) is regarded as 
string theory gauge fixing and physical quantities should not depend on it. We can think of our integration 
by parts in (|2.8p as equivalent to adding a surface term to the Wess-Zumino term in order to restore this 
gauge invariance. 

2 A monopole bundle has quantized flux. Here the number of quanta is very large in the strong coupling 
limit 2-kud = *J\f, so that it is to a good approximation a continuously variable parameter. 
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Then (|2.7[) and (|2.8|l are invariant and ()2.12p transforms covariantly under our definition 
of C. We have established that the mathematical problem of finding the D7 embedding 
has the discrete symmetries which will appear as C and P for the 2+1-dimensional 
Fermions. In the case of C this is clear, the worldsheet gauge field is dual to a 
conserved U(l) current j a = xjjj a ip for the Fermions and A — > —A corresponds to 
j a — > — j a , which is implemented by the usual C transformation of fermi fields. For 
parity, this is also the case, the Fermion kinetic term is covariant under x — > —x. 
However a 2-component Fermion mass operator iptp changes sign under this parity 
transformation, so parity takes the place of chiral symmetry in 2+1-dimensions in that 
it protects the masslessness of Fermions. Indeed, the geometric argument of reference 
|33| explains that the deviation of the angle tp from its symmetric value j corresponds 
to a separation of the D3 and D7 branes and a mass for the D3-D7 strings. At the same 
time, this deviation would violate parity. 

Our Ansatz for a solution will describe a D7-brane covering the whole range of 
(t, x, y, r, 9, <p, 9, (p) and it is for the most part determined by symmetry, 



t = a ,x =ai,y = a 2 ,z 



z(a 3 ),r = cr 3 , 6> = ai,4>= a 5 ,9 = a e ,(p = a 7 ,ip = ip(a 3 ) 

(2.15) 



we will denote the coordinate 03 by r. The two unknown functions in this embedding 
are then z(r) and ip(r). Our Ansatz for the world- volume gauge fields is 



L 



L 2 f 



2na' 



-a'(r)dr A dt + - — -bdx A dy + - — ~ ( d cos 9 A dip + d cos 9 A dip) (2.16) 

Z7TQ: Onrrv' \ I 



2na' 2 



where 6 is a constant magnetic field which is proportional to a constant magnetic field 
in the dual field theory and a(r) is the temporal world- volume gauge field which must 
be non-zero in order to have a uniform charge density in the field theory dual|f| 
With the Ansatz (|2.15|) and (|2.16p . the D7-brane world- volume metric is 



ds 2 



L 



'{-dt 2 +dx 2 
cos 2 ip(d9 2 + si 



+ dy 2 ) 
,2 a j ±2 



^(l + rV 2 + rV 2 )+ 



9d(p 2 ) +sin^(# 



(2.17) 



where prime denotes derivative by r. The Lagrangian is 



jj: = -\f(J 2 +4cos 4 V>)(/ 2 +4sin 4 V')(& 2 +r 4 )(l + r 2 i/>' 2 + r i z' 2 ) - a' 2 

+fr 4 z' + 2abc'(iP)iP' 

where, using (|2.3I) . 

XN , 



8tt 4 



2-1 



(2.18) 



(2.19) 



The factor of (2-7r) 2 in the numerator comes from half of the volume of the unit 2-spheres 
(the other factors of 2 are still in the action). Since nothing depends on (x,y,t), the 
integral over these coordinates produces the volume factor V2+1 which appears in Af. 
The Wess-Zumino term gives a source for z(r), so that, as long as the flux / is nonzero, 
z(r) will be non-zero and r-dependent. Note that, since V2+1 has dimension of length 



3 b and q are related to the physical magnetic field and charge density as b = ^j^B, q 



^■^■pso that 

AiV Z7T ' 



q 

— = TIV 

b 



1 2-Kp 
N B 



where the dimensionless parameter v is the filling fraction of N degenerate Landau levels. 
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cubed and r has the dimension of inverse length, the integral of (|2.18[) over r will be 
dimensionless, as it should be. 

Now, we must solve the equations of motion for the functions ip(r), a(r) and z(r) 
which result from the Lagrangian (j2.18[) and the variational principle. Since the La- 
grangian depends only on their derivatives and not on the variables a(r) and z(r) 
themselves, a(r) and z(r) are cyclic variables and they can be eliminated using their 
equations of motion, 



d SS 



d SS 



J(f 2 +4 cos 4 ?/-) (f 2 +4 sin 4 V>)0 2 +r 4 )rV 

x = 1 ======== /V = p z (2.20) 

dr 5z'{r) v/l + rV 5 + r 4 z' 2 - a' 2 

J (/2 + 4 cos 4 ip) (f 2 + 4 sin 4 ip) (b 2 + r 4 )a' 

= -> -i 2bc=-q (2.21) 

dr 5a' (r) ^1 + r 2 V>' + r 4 z' 2 - a' 2 



where p z and q are constants of integration. p z can be interpreted as being proportional 
to the pressure in the z-direction and, by translation invariance, for a single brane, we 
would expect that it would be zero if the brane is not accelerating. It is also clear 
that, if the brane is to reach the Poincare horizon at r — > 0, equation (|2.20p will 
make sense there only if we set p z = 0. Also, to get (|2.2ip we have integrated 2bc' to 
get the term 2bc(ijj) in the equation. We have taken the integration constant so that 
c(? — ip) = —c(ip), so that it has the correct transformation property under P, 

c(ip) = ip- f - ±sin4V> (2.22) 

The overall constant of integration in (|2.22p is, of course, a superstring gauge choice. 
Other choices would give equivalent results, but would alter what we mean by total 
charge density and would complicate the parity transformation law. Here, for simplicity, 
we make the choice given in (|2.22[) . The integration constant q is then proportional to 
the total charge density in the field theory dual. Moreover, we can solve for z' and a', 



^(P + 4 cos 4 V>)(/ 2 + 4 sin 4 ^){b 2 + r 4 ) - / 4 r 4 + {q - 26c) 2 
(26c- 9 )v / l + r 2 ^ 2 

a = - 

sj{p +4cos 4 ^)(/ 2 +4sin 4 ?A)(6 2 + r 4 ) - / 4 r 4 + (q - 26c) 2 
We must then use the Legendre transformation 

K = S- a'{r)—-r- / z'{r)- 



(2.23) 
(2.24) 



da'(r) J dz'(r) 
to eliminate z' and a'. We obtain the Routhian 

TZ = Af J dr^J(p +4cos 4 ?/>)(/ 2 +4 sin 4 ^)(6 2 + r 4 ) - / 4 r 4 + (q - 26c) 2 ^1 + r 2 ijj' 

(2.25) 

which must now be used to find an equation of motion for ip(r). Applying the Euler- 



ip' 2 



4 The parameter p z becomes important when there is more than one brane, as the branes interact with 
each other and a pressure is required to hold them together, or apart, depending on their relative orientations 
|31) . With a single brane, because of translation invariance, this pressure must vanish. 
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Lagrange equation to the Routhian (|2.25p yields the equation of motion 



ip 



1 + ip 



2r 4 [(/ 2 + 4 cos 4 ip)(f 2 + 4 sin 4 V) - / 4 ] 



if 2 + 4 cos 4 ip)(f 2 + 4 sin 4 ^){b 2 + r 4 ) - / 4 r 4 + (q- 2bcf 



4 sin 2^ [[f 2 (b 2 + r 4 ) - r 4 sin 2 2ip] cos 2^ + &gsin2?/> - 2b 2 (ip - it/ 4) sm2ip] 

(/ 2 + 4 cos 4 ip)(f 2 + 4 sin 4 V)(& 2 + r 4 ) - / 4 r 4 + (q - 2fcc) 2 ~ 

(2.26) 

where the overdot is the logarithmic derivative ip = r4-ip. In the next few Sections, we 
will discuss some of the solutions of this equation. 

3 Conformal field theory 

Let us first review the relevant solutions for the D7-brane geometry when the external 
magnetic field and the charge density vanish: b = 0, q = 0. These are described in 
references [5] and [53]- [2H|- The equation for ip(r) is 

^ 3 ; 2sin4^(/ 2 -sin 2 2^) = 

1 + V> (/ 2 +4cos 4 V)(/ 2 + 4shr» - / 4 ' 

In the large r regime, we require that the angle approaches the parity symmetric so- 
lution, ip — »• f . We also note that the constant angle ip = \ is an exact solution of 
equation (|3.27p . Linearizing about that solution yields the differential equation 



dr 



which is solved by 
where 



d \ ^, ^, §(W 2 ) 



W " f ) + 3 ^ - f ) + (/2 + 1)2 _ /4 (^ ~ f ) = ( 3 - 28 ) 




9 + i)2 _' ) y4 (3.30) 

The argument of the square root is non-negative and the exponents are real numbers if 
f 2 is large enough, 

•'' ' < 331 > 

This is the parameter regime where mass term in the linearized equation (|3.29p does 
not violate the Breitenholder-Freedman bound. In the regime || < f 2 < 1, both A + 
and A_ are positive, so that both solutions in (|3.29p go to zero at r —> oo. Moreover, 
when they are both positive, both solutions diverge at r — > 0. Thus we see that there 
is no small deviation from the solution ip = j which is itself a solution. For this reason 
we call ip = ~ an isolated solution. 

Let us discuss this in a little more detail. Regardless of its behavior at r > 0, if 
ip(r) is to have a "normalizable mode" which remains finite at r — > 0, for small values 
of r, the solution must also go to a zero of the last term in (|3.27p . When || < f 2 < 1, 
there are three zeros of the last term in (|3.27p . ip = 0, j, \ arcsin/. The exponents of 
the linearized equation in the region r — > 0, in each of the three cases are 

^ = +^ L + ^ + --- ' 7 -=l, 7-=2 (3.32) 



7T K\ K2 3 1 / 32(1 -f 2 ) 

i, = l + — + — + - ' 7±= 2 ± 2\/ 9 ' (/ 2 + l) 2 -/ 2 (3 - 33) 
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The first solution ip = would be acceptable only if both m and K2 are zero. The 
solution would then necessarily be a constant and violates the boundary condition at 
large r, where ip should go to j. If the boundary condition at infinity were different, 
so that the constant ip = everywhere were a solution, this solution would describe a 
D5-brane with flux in this geometry |42| . This is an interesting possibility, but not the 
one that we need here. 

The second solution (|3.33|) is the isolated ip = f solution. It is constant and isolated 
for the reasons that we have discussed above and it is the only solution that does not 
violate parity. 

The third small r solution is ip(r — > 0) — > ^arcsin/ which is not isolated - since 
7_ is negative and 7+ is positive, K2 is allowed to be nonzero and we need Ki = 
in order to have good behavior at small r. At large r, both constants c and m in 
equation (|3.29p are allowed to be nonzero, so the asymptotic behavior of the solutions 
is characterized by three nonzero constants, k%, m and c. Generically, for a fixed value 
of one of the constants, for example, m, both c and K2 must be tuned to obtain n\ — 0. 
Then, c and K2 are functions of m. This means that we have a 1-parameter family 
of solutions parameterized by m. This parameter is the holographic dual of a parity 
violating Fermion mass term in the quantum field theory. The ip = f solution has a 
relevant parity violating operator with conformal dimension 




32(1 - p) 

(/ 2 + 1)2 _ JA 



which can be introduced and which has coupling constant dual to the parameter m and 
expectation value dual to the parameter c|j Turning on this operator preserves charge 
conjugation symmetry but breaks parity symmetry and it corresponds to making m 
nonzero. Once m is nonzero, no matter how small, the solution will evolve to the 
same value, ip = ^arcsin/ at r = 0. This can be significantly different from the 
ip = ^ of the parity invariant solution. This has the interpretation of a renormalization 
group flow driven by a relevant operator. Note that, when / = 1, the exponents are 
(A_,A+) = (0,3). In the dual field theory this is interpreted as having the Fermion 
mass operator with dimension 3, that is, exactly marginal. In this case the angle goes 
to \ arcsin(l) — > f at small r. A numerical solution of (|3.27p us easy to find. Plots of 
numerical solutions exhibiting this behavior were given in reference [5]. 



4 Conformal field theory with charge density 

We note from the Lagrangian (|2.18p that it contains a cubic term abc'(ip)ip' ■ Remember 
that a is C and CP violating, b is P and C violating and ip' is P and CP violating. This 
implies that turning on any two of these fields will violate all three of the symmetries P, 
C and CP and, in the equation of motion, turning on two will produce a source for the 
third. In the following sections we will discuss the situation when all three are turned 
on. Before that, let us examine what happens when only one of the three is turned on. 
In the above, we have already discussed one of the cases, the situation when a and b are 
turned off but ip could be a function of r. That solution violated parity but preserved 
charge conjugation invariance. In this section, we will keep ip a constant, set b to zero 
and turn on q. This solution violates charge conjugation invariance, in that there is 
a fixed non-zero charge density, but it is invariant under parity. What we obtain is a 
sector of the conformal field theory where the charge density is fixed to a particular 

In an alternative quantization, if A+ and A_ are both greater than |, the unitarity bound for conformal 
dimensions of operators in d=2+l, the roles m and c can be interchanged [43] . 
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value proportional to q. The equation for tj)(r) is gotten from (|2.26p by setting 6 = 0, 



+ 



2[r 4 (/ 2 +4 cos 4 ^){f 2 +4 sin 4 tp) - f 4 r 4 } 
r 4 (f 2 +4 cos 4 ^)(/ 2 +4sin 4 V>) - f 4 r 4 + q 2 
2r 4 sin4^(/ 2 - sin 2 2ij)) 



= 



(4.35) 



r 4 (p + 4 cos 4 V)(/ 2 + 4 sin 4 ^) - / 4 r 4 + q 2 

We see that ip = j is a solution of this equation. The remaining embedding functions 
can be found from by integrating the expressions in (|2.23|) and (|2.24p with 6 = and 

f 2 



dr 



o 



x /(2/ 2 + l)f4 + g 2 

-5 



(if — == — 

v/(2/ 2 + + q 



(4.36) 
(4.37) 



Note that the range of these function is finite. 

The chemical potential is given by difference in the value of the gauge field at r 
and r = (where we remember the normalization in (|2.16[) ) 



u = — — a(oo) - a(0) = -— r O -, - srgn(g) 

' 2tt ; WJ 2tt 4(2/2 + 1)3 L 4 4 J 6 v ; 



(4.38) 



B[s,t] = r[s]r[t]/r[s + t] is Euler's beta-function. Now, remembering that the charge 
density is defined by p = q, we can write the expression for the charge density as 
a function of the chemical potential, 



p = -sign(/ti) 



16^2/ 2 + 1 N 
■kB[\,\] 2 V\ ' 



(4.39) 



For a generic value of f 2 in the range of interest, say f 2 = 1/2, 
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p sign(|U) —=Np 2 with / 2 = 1/2 

v A 

One might compare this with the free field theory where 

Po = - p£pr J d 2 kd(p - \k\) = p 2 ~ -.080iV M 2 

The scaling with /i 2 is a consequence of conformal invariance. In a theory where the 
U(l) current is conserved, the Debye screening mass (the inverse of the Debye screening 
length) can be derived from the charge density-chemical potential relationship by taking 
a derivative of the charge density by the chemical potential |44j. 



(4.40) 



(4.41) 



dp 



32y/2/ 2 + l TV _ 8 (l + 2/ 2 )i y/N 



1^ 



(4.42) 



It is interesting that coefficients of Np 2 in the strongly coupled theory and the free 
field theory and the Debye screening lengths differ by a factor of A= which can be 
significant in the large A limit. The Debye screening length as a function of gate voltage 
is a quantity which could in principle be measured in a relativistic condensed matter 
system such as graphene. The coupling constant could also be varied, in principle, by 
changing the dielectric constant of the vicinity of the material. Such a measurement 
would be an interesting test of the hypothesis that the relativistic material could be in 
a strongly correlated state that is described by a conformal field theory similar to the 
one which solves the D3-D7 model. 



10 



5 Conformal field theory in a magnetic field 

The equation for the angle ip with 6^0 and q = is 



i + i> 2 



2r 4 [(/ 2 + 4 cos 4 ip) (/ 2 + 4 sin 4 V) - f 4 



{f 2 +4 cos 4 ^)(/ 2 +4sin 4 V>)(6 2 +r 4 ) - / 4 r 4 + (26c) 5 



4sin2?/> [[/ 2 (6 2 + r 4 ) - r 4 sin 2 2ip] cos 20 - 26 2 (0 - tt/4) sin 20] 
(/ 2 +4cos 4 V')(/ 2 +4sin 4 ?A)(6 2 +r 4 ) - / 4 r 4 + (26c) 2 



(5.43) 



This equation is solved by tp = ^. The magnetization as a function of field is given by 
the derivative of the vacuum energy by the field, 



M 



d_ 

dB 



F 



1 



8 



V 2 +i dB 



1Z 



(5.44) 



where we have recognized that the energy as a function of field is simply the negative 
on-shell Routhian, which we can find by setting q = 0, ip' = in (|2.25p . We obtain 



K=N drV(.f 2 + l) 2 (6 2 +r 4 )-/V 



(5.45) 



The above integral diverges. It can be defined by subtracting the 6 = energy from 
it. Then the result is finite. Then, the derivative by the field renders the remaining 
integral finite, 



M 



\N ( 2tt 

8^ v7x 



B, ir - + 

o V(/ 2 + l) 2 (l + H)-/M 



sign(S) 



M = -■ 



(/ 2 + l)UJiV 
2(27r)l(2/ 2 + l)i 



■B 



i i 



B sign(B) 



M (0.06) X^nVb sign(B) (with f 2 = 1/2) 



(5.46) 



The sign, M ~ — sign(i?) is a result of diamagnetism, which is what is expected for 
electrons in a magnetic field. We can compare this with the result for free Fermions. 
There, all negative energy states are occupied and their spectrum in a magnetic field 
is that of relativistic Landau levels, given by E n = -y/2\B\n, with n = 0, 1, 2, .... The 
degeneracy of each Landau level is The ground state energy is 



27T 



\B\iN 
V2tt 



C(-l/2) 



where we have defined the infinite sum by zeta-function regularization and the value of 
the zeta-function is ^(—1/2) = —0.21. The vacuum energy is positive and 



M = -- 



Sir 2 



(0.21)Ny/\B\ sign(B) - -0.07NVB sign(S) 



(5.47) 



It is again interesting that the diamagnetism of strongly and weakly interacting Fermions 
differ by a factor of A*. 

Now, we observe that, with a finite magnetic field, at r — > 0, the angle "0(0) must 
solve the equation 

f 2 cot 2^(0) - 2(0(0) - tt/4) = (5.48) 

The only solution of this equation is f/>(0) = ?, which fixes the value that 0(r) must 
take at r — s- 0. However, if we study the linearized equation in this region, we find 
that, unlike the case that we studied in the previous sections where the constant i(j = \ 
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was an isolated solution, here, it is not isolated. There are solutions close by where 
ip(r) depends on r and approaches j when r goes to zero and infinity. To see this, the 
linearized equation in the small r region is 



and it has solutions 



7T l + I /l i 32 ,, I-i [\+~- 

ij}{r) = - + Kl r 2 + 2 V +n 2 r 2 2 V +T ^+... (5.50) 

We see that the solution (|5.50p converges to ^ if we choose the constant k 2 = 0. This 
requires tuning the asymptotic behavior of the solution at large r, which we shall exam- 
ine in more detail in the following. In this way, we see that, when there is a magnetic 
field present, there are solutions of the equation for ip(r) which are infinitesimally close 
to the constant ip = ^ for all values of r. This is in contrast to the case without a 
magnetic field where any infinitesimal deviation of ip(r) from j at large r led to a non- 
infinitesimal deviation at small r. In the present case with magnetic field, we can study 
solutions which deviate from j by a small amount perturbatively. We shall spend the 
remainder of this section studying the properties of these solutions. 

Before we begin, we note that, as soon as ip( r ) deviates from the constant j, the 
solution violates all three of the discrete symmetries, parity, charge conjugation and 
CP. That it violates C and CP can be seen from the equation (|2.24p which requires a(r) 
to be non-zero, at least for intermediate values of r. However, we are still free to set 
the parameter q, and therefore the charge density to zero. In the next Section, we will 
examine what happens when we turn on an a q of infinitesimal magnitude. 

We begin by considering a solution which deviates from the constant by a small 
amount, 

V>W = J + 0(r) (5.51) 
and Taylor expand (|2.26p to first order in </> and q, we obtain 

[b 2 (l + f) 2 + r 4 (l + 2/ 2 )]rV + 2[6 2 (1 + ff + 2r 4 (l + 2/ 2 )]r^'+ 

+8[r 4 (l - f) - b\l + f 2 )}4> = -4bq (5.52) 

which we shall now try to solve over the full range of r. Note that (|5.52p reduces to our 
small r linearization in (|5.50p at r — > 0. The transformation 



pF \ 1 + 2f 2 r 4 

y( x ) = y^r^( r ) i x = y 1+ (i + /2) 2 ft2 ( 5 - 53 ) 

puts (|5.52p in the following form, which in the case of q = is the standard form of the 
associated Legendre equation 

(l- a V-W+ + 1) - t^Zs) y - -I" ow 1 , . ^1 (^i) 8 (5-54) 

where the parameters are 



1-X 2 J 6(l-a; 2 )(l + / 2 )t VI + 2/ 




50/2 ' 23 u=^f+l_ (5 . 55) 



16(1 + 2/2) ' 4 V /TT7 2 ' 

For future reference, we note that u is always positive but v can change sign. One can 
write v in terms of A± as 




, I-/ 2 
'2/2 + 




(5.56) 
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4 2 

The Legendre equation has the general solution 



1 



ClP 



(l + 2/ 2 ) r 4 
(1 + J 2 ) 2 6 2 



c-iQ 



u,u, Wl + 



(1 + 2/ 2 ) 
(1 + ,/ 2 ) 2 b 2 



(5.57) 



(5.58) 



The associated Legendre functions are given by Hypergeometric functions as below 



ri-u 



1 + X 



-1 + :r 



2F1 



—v, v + 1; 1 — u\ 



1 - x 



(5.59) 



r[;/ + u + 1] (-1 + x 2 ) 

2 V + 1 T[v + 'i/2] 



2\u/2 



-%r 1 



u + i/+lu + i/ + 2 n . 1 

— 5 — > — 9 — ^ + 3 /2;-2 



(5.60) 

in which we have defined phases so that the functions are real for x > 1. Both of 
the solutions (|5.59p and (|5.60p diverge at r — > (which corresponds to .t — > 1 + ) and 
converge to zero at r — > 00 (which corresponds to x — > 00). 

To find a solution which is finite at r — > 0, we need to adjust the two constants, C\ 
and C2 in (|5.58p so that the singularity cancels. Since r — > corresponds to x — > 1 + , 
we shall need to study 2Fi[a, b, c, z] in the regions where z = for P and z = 1 for Q, 
respectively. It is easy to study P near x —> l + . From the definition of Hypergeometric 
functions. 



2F1 [a,b,c,z] = ^2 

n=0 



(a) n {b) n z* 



(a)„ = 



T\a + n] 



(c)„ n\ ' v "" r[a] 
we see that the leading term in the series is normalized to one, 

2 F 1 [a,b,c,0] = 1 

Then, from (|5.59[) , we see that the asymptotic behavior of P is given by 



1 



r[i - u] 



-l + x 



x ~ r 



(5.61) 



(5.62) 



(5.63) 



Now, let us consider Q. We note that, for the values of u and v if interest, Q 
is divergent at z = 1. Since the asymptotic small r behavior of the solutions of the 
Legendre equations exhibit only one divergent solution, we know from the outset that 
the divergent nature of Q must be similar to that of P. 

Euler's formula is an integral representation of F which is valid when 5R(c) > and 
R(6) > 0, 



2 Fi[a,b, c, z] 



Using it, one can show that 



r ^ r r[6] / t b - 1 {t-iy- b - 1 {i-tz)- a dt 



T[c-b] 



2 Fi[a,b,c,z] = (1 — z) 



c—a—b 



F\ [c — a, c — 6, c, z\ 



Using the same representation one can show that if 5ft(c — a — b) > 

T-1 r t n T[c]T[c-a - b] 

2 Fl fl, b, C, 1 = — : r— 

1 [c — ajl [c — b\ 

for $t(b) < the above result can be generalized using analytical continuation. 



(5.64) 



(5.65) 



(5.66) 
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In our case of interest, Q in Eq. (|5.60p has 

, u + v + 2 

b = > 



— a — b = —u < 



(5.67) 
(5.68) 



which is not in the domain of applicability of (|5.66[) . If we use the transformation (|5.65[) . 
we can change the arguments of F so that (|5.66p can be applied. Then we get 



where, now 



The result is 



2-Fl [a, 6, c, z] = (1 - z) u 2 Fi [c - a, c - 6, c, z] 

c— (c — 6) — (c — a) = a + b — c 
= u > 



(5.69) 

(5.70) 
(5.71) 

2 „+i \2{-i + x)J r[2±f±i]r[2±|±2] ' (5 - 72) 

which shows the expected divergent behavior. Then the regular solution turns out to 
be (|5.58p (which we recopy here) 



2 Fi[a, b, c, z\ ~ (1 - z 
and we find the asymptotic behavior of Q, 



_ M r[c]r[a + b-c] 
r[o]r[6] 

r[«] 



/2 



1 



ClP 



V, u, \ / 1 + 



1 + 2/ 2 r 4 



c 2 Q 



(1 + py 62 
in which c\ and c 2 are related by 

c 2 _ 2^+ u+1 r [ !i± ^ ±1 ] r f 34 ^^] 

Cl 



1 + 2/2 r 4 

(1+7^62" 



(5.73) 



r[u]r[i/ + u + i]r[i - u] 

We have now found an acceptable non-singular solution in the small r — > regime. 
It is given by equation f|5 . 58[) with the constants related by (|5.73p . We now need to find 
what this implies for the large r behavior of the solution. In particular, it should fix 
the relationship between the coefficients of the two asymptotic behaviors of the solution 
at large r which are displayed in equation (|3.29|) . that is, it should fix the relationship 
between c and m. For large r one gets 



x = 1 



1 + 2/2 r 4 ^1 + 2/2 r 2 



^ (l + /2) 2 62 (1 + /2) b 

Using the same integral representation (|5.64p . one can show that 

z 



2 F 1 [a,b,c,z] = {l-zy a 2 F 1 
For P in the regime where - ~ we get 



a, c — b, c, 



z - 1 



(5.74) 



(5.75) 



1 



1 



r[i - u] \ 

1 fx 



—v, v + 1; 1 — u) 



l-x 



T[l-u] V2 
1 fx 



1 











-) 




— ) 2Fl 


—v, —u 


X I 




xJ 





(5.76) 



r[i -u}\2 



V + u 



X 



2 d 



2 i*i [—v, —u — v, 1 — u, 1] — 9F1 [—v, —u — v, 1 — m, x] 

x dx 



(5.77) 
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Because 1 > v > — |, then we find that 1 — u — (— u — u) — (—v) > 0, which allows us 
to use (|5.66[) . Using the identity (|5.66|) . we find that 



2-Fi [—v, —u — v, 1 — it, 1] 



r[i - u]r[i + 2v\ 



T[l-u + v\T[l + v] 
Another identity that we shall use is the following identity 

2 F 1 [a,b,c,z] = — 2 Fi[a+l,b + l,c+l,z] 

az c 

Using this identity, we see that the second term in (|5.77p can be written as 



— 2 Fi \—v, —u — v.l — u, x] 
dx 



v(y + u) 
l-u 



2-F1 [1 — v, 1 — u — v, 2 — it, x] 



(5.78) 



(5.79) 



(5.80) 



Now, we should be more careful because for > v > — i, we can not use the integral 
representation (|5.64p . In this case c — a — b = 2v < and the (|5.66[) can not be used. 
We shall therefore consider the two cases v > and v < separately: 
For v > 0, using (|5.80p and (|5.66|) . we find out that 



dx 

Then one has 



Fi [—v, —u — v, 1 — u, x] 



v{v + u) V[2-u]T[2v] 
l-u r[l-u + z/]r[l + i/] 



(5.81) 



in which the terms that fall off as x v ~ Y get canceled by each other (as one expects from 
the asymptotic solutions of (|5.54[) for large values of x). 

To deal with the v < case we can not use the Taylor expansion directly, as the 
derivative of F at x = 1 diverges. From the asymptotic solutions of (|5.54p . we know 
that to the two highest orders, P can either go like x" or x~ 1 ~ u . The x v part is easy to 
find, we can just use the (|5.76[) . and look at x — > oo. If we subtract the x v term from 
P, because x~ l ~ v vanishes at infinity, the result would vanish at infinity. Considering 
this fact, we find the following limit 



lim 

x— >-+oo 



iFi [—v, -u - v, 1 - u, 1 - -] - 2-Fi [—v, —u - v, 1 - u, 1] 



-2i/-l 



d tF\\—v, — u—v,\ — u,\- 



lim 



dx 



(5.83) 



dx 



which can be used to find the coefficient of , 



term in P (the extra power of v 



is because of presence of x v in (|5.77[) ). Using (|5.80p . (|5.65p and (|5.66p . we find that 



d F 

ax 



-v, —u — v, 1 — u, 1 

x 



v{v + u) r[2 - u]T[-2v] 
l-u r[l - v]T[l -u-v\ \x 



(5.84) 



in which comes from the chain rule used to change the variable of differentiation. 
So we find that for v < 0, 2-Pi [ — v, —u — v, 1 — u, 1 — — 1 can be written as 



iF, 



-v, —u — v, 1 — u, 1 



2-F1 [— i/, — u — i/, 1 — u, 1]- 



+ u) 



r[2 - u]r[-2v] 



(1 - u){2v + 1) T[l - v]Y[l -u-v) 



2u+l 



(5.85) 
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Plugging in the above results in (|5.76[) . we find the following asymptotic expansion 
of P 



T[l + 2u] 



r[i-u + i/]r[i + i/] 



(§)" 



1/ + u 



v+1 



v(y+u)- 



Y[-2v-l] 



T[l - v]T[l -u-v] 
(5.86) 

where we have used the fact that T[x + 1] = xr[x]. 

We neglect the term, as it belongs to next orders. We then have the following 
simplified equation of P for v < 



T[l - u + v]T[l + u] V2 



v(y + u) 



T[-2u- 1] 



r[l - u]T[l -u-v] 



(5.87) 



Around v = — ^ the normalizable/un-normalizable nature of the two terms switches. 
Later, we will explore the behaviour of the solution around v = — A, and see if there is 
anything special about this point. 

Now we shall study Q near the boundary, r — > oo. We Taylor expand Q by using 
the fact that 2-F1 [a, b, c, 0] = 1 and (|5.80|) . which result in 



u + iv + l u + v + 2 . 1 



1 



+=0 



1 (u + i/ + l)(u + t/ + 2) 
x 2 4(i/ + 3/2) 



(5.88) 



we find that 

0rT[i/ + « + l] 1 



2"+ 1 r[i/ + 3/2] 2^ +1 
yiT[i/ + it + l] 1 
2"+ 1 I> + 3/2] x^ 1 



u 

2x2 
u 

~ 2tf 



1] 2*1 



m + z/+1u + z/ + 2 o . 1 
-51/ + 3/2; -j 



2 2 
1 (u + i/+ l)(u + t/ + 2) 
z 2 4(i/ + 3/2) 



(5.89) 



We expect to hnd the asymptotic behavior (|3.29[) . that is, <p(r —> 00) 
Here, the correct powers of r are found by 



1 , / 50f 2 -23 



4(1+2/*) ^ r A_ 



(5.90) 
(5.91) 



in which the extra factor of r comes from definition of <p in (|5.58D . 

The first asymptotic term (with m), which decays slower, always comes from P. 



Y[l + 2v] 



T[\-u + v]T[\ + v\ \2 



iy ( ^T+2p \ 
1 1 + /2 b 



Cl 



(5.92) 



The term is coming from the fact that 4> is defined in terms of r but the 

Legendre Functions are in terms of x. 

Depending on the sign of v, the second asymptotic term (with c) may come from 
either Q or from both of P and Q. If v > 0, the second term in P goes at most like 
x"~ 2 and first term in Q goes like x~ v ~ x , because v < 4, the term coming from Q is 



the dominant one. (We recall that A + = 2v + | and A_ 
Looking for the coefficients, one finds that 



21/.) 



rp±|±i]r[«± i ±2] / 1 + /2 



r[u]r[i- u ]r [!/ + §] VVT+2? 



Cl 



(5.93) 
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Table 1: Special values of f 2 



To 

/ 2 






Zo 

50 


O O 

2' 2 


stability bound 


1 

2 


2,1 


classical dimensions 


9 
14 


5 1 
2' 2 


unitarity bound 


0.74 


m = 0, c t^O 


chiral symmetry breaking 


1 


3,0 


marginal 



Table 2: Table of special values of f 2 . For stability of the system to fluctuations, it is 
necessary that f 2 > |l. When / 2 < yj, both dimensions satisfy a unitarity bound and there 
exists an alternative quantization where the mass and the condensate can be interchanged. 
At f 2 ~ .74, magnetic catalysis of chiral symmetry breaking is can take place. When f 2 = 1, 
the mass operator is exactly marginal. 



then 



_£ = _ 2 „+„ f i + / 2 A 2 ' +1 r[i- u + „] r[i + ,]r[HMi]r[^] 

m ^ v/i + 2/ 2 ; r [i + 2v] i>] r[i - u] r [i/ + 1] 1 ■ j 

If ^ < 0, we use the expansion we found for P and Q at a; = oo, which results in 



2" 



r[i - 2*y] 



(v + it) 



r[i - v\ r[i 



u] 2^ + 1 



ci 



r [h±|±i] r [gfg] 



r[u] r[i - u] r [i 



(5.95) 



then 



( 



l+f 



t) 2 



2" 



rp-gt/i 



rri-firfi-v-u] 2i/+i 



r[«]r[i-«]r h/+j 



r[l+2;y] 



r[i-u+Hr[i+H 



(5.96) 



This change of the ratio — which occurs when the flux f is varied so that v crosses 
zero resembles a quantum phase transition. One should note that none of our restrictions 
about parameters in Hypergeometric functions were violated. Use of the linearized 
equation is justified here a postiori as the solutions remains small. 

We note one other interesting behavior which occurs as / is varied. When it achieves 
a value f 2 ~ .74, so that 1 — u + v = 0, the ratio in (|5.96p seems to diverge. This 
behavior is depicted in figure . We interpret it as the condensate c remaining non-zero 
as m is put to zero, that is, spontaneous symmetry breaking, the presence of a nonzero 
condensate when the source is put to zero. In addition, at this point, the sign of the 
condensate versus the sign of the magnetic field flips, for f 2 > .74 the sign of c and the 
sign of m are identical, for f 2 < .74 the signs are opposite. 

We have found that the conformal dimensions of the operator which is dual to 
fluctuations of the angle tp lie in an interesting range when |S < f 2 < 1. A summary 
of some special values of f 2 in this range are given in Table [5] 
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Figure 2: The ratio of condensate to mass parameter ^ is plotted versus the internal flux parameter f 2 . 
The plot exhibits a discontinuity at v = and a singularity at 1 — u + v = 0, the value where the 
condensate changes sign, and also where spontaneous breaking of chiral symmetry can take 
place. 

6 Finite density in presence of magnetic field 

As we have remarked in the previous sections, the equation of motion (|2.26p has the 
constant solution tp = it /A only when either b = 0, q = or both. When both 6 and 
q are nonzero, there is no constant solution. When q and b were set to zero, we saw 
that the constant solution tp = ? was isolated. As soon as we turn on a non-constant 
behavior of ip(r) at r — > oo, the solution near r — > deviates by a large amount from 
j. Then, in the previous section, we found that when b is nonzero, but we still tune q 
to zero, the constant solution is no longer isolated and we explicitly found a solution 
which differs from the constant by a small amount for all values of r. In this Section, 
we shall examine what happens when we turn on a small value of q. 

In the small r region, the non-constant solution must go to a zero of the last term 
in (|2.26p . That is, it must solve the equation 

f cot 2^(0) + | - 2(>(0) - tt/4) = (6.97) 

When | = 0, this equation is solved by "0(0) = f, and the constant solution is allowed 
and, as we saw in the previous section, even when it is not constant, ip{r) approached 
j when r approaches either or oo. Now, when | is nonzero, this is no longer the case, 
there is no constant solution and the non-constant must approach a zero of (|6.97p which 
is no longer j . If we consider the case where | is small, it is solved by 

HO) = (6.98) 
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It is consistent that the deviation of tp(r) from j is of order |. If we consider 

(r) - \ + </>(r) (6.99) 

and Taylor expand (|2.26p to first order in 4> and q, we obtain equation (|5.52p . which we 
recopy here for the reader's convenience: 

[b 2 {l + ff + r 4 (l + 2/ 2 )]r 2 0" + 2[6 2 (1 + / 2 ) 2 + 2r 4 (l + 2/ 2 )]r<//+ 

+8[r 4 (l - f) - 6 2 (1 + f 2 )]4> = -4bq (6.100) 

The solution of the (|6.100|) is given by 

/>oo 

4> = 4>h + dyG(x-y)f(y) (6.101) 
Jo 

in which <pH is the general homogenous solution and G(x — y) is the Green function. 
The Green function is given by 

G(x, y) = g{x)h{y)Q(y - x) + g(y)h(x)Q(x - y) (6.102) 

in which g(x) is the regular solution at origin, h(x) is the regular solution at the bound- 
ary and is the step function. Having only one regular solution at the origin, we are 
forced to choose the linear combination of P and Q that we found in (|5.73l) . Both of 
P and Q are regular at the boundary and vanish. We specify the solution by specifying 
the boundary condition at r — oo so that the Green function part does not contribute 
to the coefficient of term that dies slower, to; so, we choose h(x) to be Q(x). We are 
just left to find the right normalization of our Green function so that 



(1 - x 2 )G(x, y)" - 2xG(x, y)' + [v(v + 1) - j^-^j G(x, y) = S(x - y) (6.103) 
in which the primes are differentiation with respect to x. This equation can be written 



^((1 - x 2 )G(x, y)') + (v(v + 1) - G(x, y) = S(x - y) (6.104) 

Integrating the above equation, we get the following constraint 

(1 - x 2 )G(x, y)'\ v+t - (1 - x 2 )G(x, y) V, = 1 (6-105) 
which can be written as 

(1 - x 2 )W[g(x), h(x)} = 1 (6.106) 

in which W[h(x), g(x)} is the Wronskian. The (|6.104[) has a form of Sturm-Liouville 
equation. It is a well-known property of such an equation that the Wronskian of the 
solutions is given by 

W[h(x),g(x)] ex 1 (6.107) 
(1 — x l ) 

which ensures that we can satisfy the constraint in (|6.106|) . We just need to evaluate 
the Wronskian at a given point to fix the normalization. Substituting the asymptotic 
expansions of P and Q, we get 

win pi - r^+V^ + n r(i + 2i/) i> + M + i) 0F 

W[Q,P]-(-) (2u + l) T[1 _ u + vW + v) r^ + 3/2) {1 _ X 2 } (6- 108 ) 
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Now we can use (|6.108|) to find the Wronskian. Using the linear combination (|5.73p as 
g(x), we find that the Wronskian is given by 



W{Q,g] = (-) 



2' r(l - u + u)T(l + v) r(i/ + 3/2) (l-x 2 ) 

= 1 c(rhj (fU09) 

So, we need to define h(x) in the following way to get the right normalization 

h(x) = CQl{x) (6.110) 
Putting all the parts together, we can finally write the complete expression of ifi 



i> = -r+<t>H ~ - h -^=Ql{x) [ dy (p( v , u,y) + ^ Q(v, u, y)) -± — T ( ^— j, 

4 b y/rb-i Jo V ci / (l-2/ 2 )(l + / 2 )" V 1 + 2/- 

(6.111) 

in which tpn is the solution to homogeneous equation and the rest are as listed below: 



/ 1 4- 2 f 2 r 4 

'i 1 + iTT?F? «" 12) 



(6.113) 



c 2 _ 2 v+u+1 r [H+^+i] r [2±|±2] 
^ _ r[u]r[z/ + u + i]r[i-u] 

L = /V+i T(2 + 2^) ViT^ + ^ + l) 

c 4 J r(i - « + u)r(i + v) r(i/ + 3/2) 1 j 

In this section, we have derived an explicit solution of equation of motion in presence 
of magnetic field and finite density. This solution shows finite deviation from constant 
solution. 



7 Conclusions 

In this Paper, we have examined a regime where a mass operator, magnetic field and 
charge density can be turned on and the modification of the D3-D7 conformal field 
theory can be studied perturbatively. In Sections 5 and 6 we solved the leading per- 
turbation explicitly. For this solution, we can compute the relationship between the 
coefficient of the mass operator and the chiral condensate for interesting values of f 2 . 
This function is depicted in figure^ It has the interesting feature that, at f 2 = .74, the 
relationship becomes singular. We interpret the singularity as a signal of spontaneous 
symmetry breaking. 

Moreover, as a byproduct, we have found holographic expressions for the Debye 
screening length as a function of fermion density and of the diamagnetic moment of 
the ground state as a function of magnetic field. They both exhibit a dependence on 
the charge density and magnetic field consistent with scale symmetry. They also have 
a mild dependence on the coupling constant which could in principle be interesting to 
compare with a real material, such as graphene where these quantities can be measured. 
It would be interesting to check the functional form in order to test the hypothesis that 
graphene is scale invariant. It would also be interesting to see whether these quantities 
in graphene differ in any significant way from the free field expressions, and whether 
the deviation is in the same direction as predicted by the strong coupling formula. 
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